The ^^0(a,7)^^Ne reaction rate and the stabihty of thermonuclear 

burning on accreting neutron stars 
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ABSTRACT 

Neutron stars in close binary star systems often accrete matter from their companion stars. 
Thermonuclear ignition of the accreted material in the atmosphere of the neutron star leads to a 
thermonuclear explosion which is observed as an X-ray burst occurring periodically between hours 
and days depending on the accretion rate. The ignition conditions are characterized by a sensitive 
interplay between the accretion rate of the fuel supply and its depletion rate by nuclear burning 
in the hot CNO cycle and the rp-process. For accretion rates close to stable burning the burst 
ignition therefore depends critically on the hot CNO breakout reaction -^^0(a, 7)^^Ne that regu- 
lates the flow between the hot CNO cycle and the rapid proton capture process. Until recently, 
the ^^0{a, 7)^^Ne reaction rate was not known experimentally and the theoretical estimates car- 
ried significant uncertainties. In this paper we perform a parameter study of the uncertainty of 
this reaction rate and determine the astrophysical consequences of the first measurement of this 
reaction rate. Our results corroborate earlier predictions and show that theoretically burning 
remains unstable up to accretion rates near the Eddington limit, in contrast to astronomical 
observations. 
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nuclear reactions — stars: neutron 



Recently, many grouDS (IWooslev et al. 


2004 


Amthor et al. 2006; Cooler & NaravanI 


2006a 


Fisker et al.l 20061 20071: Roberts et all 2006) have 



considered the nuclear reactions of the explo- 
sive thermonuclear runaway that lead to a type 
I X-ray burst on the surface of an accreting 
neutron star (for reviews, see iBildstenI Il998l: 
Strohmaver fc Bildstenl20"ol . The i50(a, 7)i9Ne- 



reaction is particularly interesting because it 
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is the gateway between the hot CNO cycles, 
which are associated with stable burning during 
the quiescent phase, and the rp-process, which 
i s associated with the t hermonuclear runaway 
()Wallace fc Wooslevlll98l[) . 

When hydrogen and helium accrete onto the 
neutron star and advect into the atmosphere, hy- 
drogen burns via the hot CNO cycle ^^C(p, 7) 
"N(p,7) ^^0{l3+ve) ^'*N(p,7) i^0(/3+i/e) "N(p, a) 
^^C The hot CNO cycle is temperature inde- 
pendent (/3-limited) and therefore its rate de- 
pends only on the abundance of ^^O and ^^O. 
Using a ti r ne-dep endent X-ray burst model, 
Fisker et al.l ( 20061 ) investigated the importance 
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of the ^^0{a,j) ^^Ne breakout reaction rate and 
showed that a lower breakout rate restricts the 
outflow from the hot CNO cycle so that the hot 
CNO cycle processes hydrogen to helium at a 
faster rate. This increases the energy output and 
thereby raises the temperature of the entire en- 
velope so that '*He burns to ^^C further out in 
the atmosphere, which in turn increases the ef- 
ficacy of the hot CNO cycle. This moves the 
hydrogen burning front (defined by Xh — 0) 
further out where ^He burns slower. The en- 
ergy generation rate therefore decreases until the 
hydrogen burning front advects down to a suffi- 
cient depth and restarts this cycle. For a global 
accreting rate of M = 1 x 10^^ gs^^, this re- 
sults in the non-bursting oscillatory luminosity 
found in time-dep endent model simulations by 
Fisker et al.l (l2006h and later in the two-zone 
model of (|Cooper fc NaravanI liooebl ). The re- 
sult of this burning behavior is a copious produc- 
tion of ^^C ( Fisker et al. 20061 ) which theoretically 
would provide the r equired fuel concentration t o 
explain superbursts (jCumming fc BildstenllioOlh . 
Unfortunately for this theory, astronomers fre- 
quen tly observe X-ray bursts at this accretion rat e 
(e.g.. lCornelisse et al.lliool ICallowav et al.ll2006[ ). 
which suggests that the aforementioned burning 
cycle does not occur. This leads to the astrophysi- 
cally based conclusion that the true ^^0(a, 7)^^Ne 
reaction rate must be higher tha n the lower limit 
calculated bv lFisker etal] (|2006l ). 

In the first part of this paper, we therefore 
perform a parameter study of the ^^0(a,7)^^Ne 
reaction rate as a function of the accretion rate 
to investigate the impact of uncertainty of the 
reaction rate. The results of this study cor- 
roborate earli er estimates of t h e crit ical accre- 
tion rate by 'Fushiki & Lamb' ("I OSTI): 



Bildsten 



il998)iFisk ei - ct al . ( 2003): Wooslev et al. (12004 ) ; 
Heger et al.l (|2005l ): iCooper fc NaravanI (|2o'o6a[ l. 



It also stren gthens the astrophysically based con- 
clusion of Fi sker et al.l ( 2006) that set a lower limit 
on the ^^0(a,7) ^^Ne breakout reaction rate 

In the seco nd part, w e use t he recently mea- 
sured rate of iTan et al ( 2007 ). which includes 
an experimentally determined lower limit and re- 
peat a nd improve the calculations of 'Fiske r et al.l 
(l2006l) . We also take advantage of the significantly 
reduced uncertainty to determine the uncertainty 
range in the transition accretion rate between sta- 



ble and unstable burning. 

We discuss the computational model in S|2] 
The parameter study is described in Sj3] which 
is followed by a discussion of the implication of 
Tan et al.l (j2007l )'s new rate in ij4]and the conclu- 
sion in |j5l 

2. Computational Model 

The impact of the ^^0(a,7)"'^^Ne rate has been 
investigated in the framework of a dynamical and 
self-consistent spherically symmetric X-ray burst 



model that has been used also in Fisker et al 



(|2005al |bl. [200il2007l ). This model couples a modi- 



fied v ersion of a general relativi stic hydrodynamics 
code (jLiebendorfer et al.ll2002l) with a generic nu- 



clear reaction network (jHix fc ThielemannI 119991) 



using the operator-split method (jHenvev et al 
19591 ). The nuclear reaction flow and the con- 



ductive, radiative, and convective heat transport 
are computed in a general relativistic spherically 
symmetric geometry. 

The reaction network comprises the same 298 
isotope network between th e valley of s t abilit y 
and the proton drip line as iFisker et al.l (|2006[ ). 
Except for the ^^0(Q;,7)-^^Ne-rate, all the proton-, 
neutron-, and alpha-induced reactions are adopted 
from the latest version o f the REACLIB li- 
brary (ISakharuk et al. II2OO6I) which was also used 
in Weinberg et al. (2006). 

The radiative, conductive, and convect ive heat 
transport is treated in the formalism of ' Thornd 
(|l977l ). The radiative opacities due to Thompson 



scattering a nd free-free absorpti on are calculated 
according to Schatz et al. I (|1999[ ) and the conduc- 
tivities for electron scattering on electrons, ions, 
phonons, and impurities are based on the work of 
Brownl (I2OOOI) . The accreted matter is assumed to 
be fully ionized on impact. We use an arbitrar- 
ily relativistic and arbitrarily degenerate equation 
of state to describe the electron gas. Due to the 
shorter quantum wavelength of the nucleons, the 
nucleon gas behaves as an ideal gas. 

The model code tracks energy transport with 
high precision and takes into account the heat 
transport betwe en the atmosphere and the neu- 
tron star core (jBrownl |2004 ). The thermal en- 
ergy released from the core is O.llMeVnuc""'^ for 
M = 1.12 X 10i'^gs"\ but it generally depends 
on the accretion rate, the composition of the X- 
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ray burst ashes, the temperature of the envelope, 
and the mass and radius of the core. The neu- 
tron star core boundary corresponds to a pressure 
of P = 7 X lO^^dyncm^^ and the core itself is 
characterized by a mass of IAMq and a radius 
of 11.06 km leading to a redshift of 1 + z = (1 — 
2GM/Rc^y^/'^ ^ 1.27 and a gravitational accel- 
eration of 5 = (1 z)GM/R^ = 1.8 X 10"cms-2. 
Similarly, the upper atmosphere is descri bed by a 
relativistically corrected grey atmosphere ( Thornd 
19771: IWeiss et al.ll2004l ) using a 4th order Runge- 
Kutta method to numerically integrate the hydro- 
static heat and pressure equations from the model 
boundary out to P = 10^* dyncm"^. 

The computation followed a sequence of accre- 
tion and burst phases until a limit cycle equilib- 
rium was reached. This was performed for dif- 
ferent ^^0(a,7)^^Ne rates for a range of accretion 
rates as described in the next section. 

3. Parameter study of ^^0(a,7)^^Ne 

Increased availability of computing power has 
significantly expanded our ability to perform pa- 
rameter studies of full one-dimensional X-ray 
burst simulations. Fig. [T] shows the results of 
a parameter study o f 72 simulations where th e 
i50(a,7)i9Ne rate of ICaughlan fc Fowled (Il988l) 



which is based on the work of iLanganke et al 
( 19861 ). has been scaled linearly with a param- 
eter, /. Similarly the accretion rate has been 
scaled linearly with a parameter, /, where 1=1 
corresponds to the Eddington limit M^^^ = 
1.12 X lO^^gs-^ The local Eddington limit is 

mgdd. = ^%dd./(4^-R') = 6.62 X lO^gcm-^s"! 
assuming that accretion is spherically symmetric. 
The simulations were computed on an 8 x 9 grid 
corresponding to the u (bursting behavior) and s 
(non-bursting behavior) letters in Fig. [TJ The grid 
covers the range between H/He ignited bursts and 
stable burning. In the following subsections we 
describe the different outcomes of this study. 
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Fig. 1. — The figure shows a log-log plot of the 
i50(a,7)i9Ne rate of ICauehlan fc Fowleij (|l988l) 
scaled with a parameter, /, and the accretion 
rate scaled linearly with a parameter, /, so M = 
/ X 1.12 x lO^^gs"^. Each letter represents a sim- 
ulation with u showing bursting behavior and s 
showing non-bursting behavior. For reference, the 
solid line shows the rate of the new measurement 
for r = 4 X lO^K. The light and dark shaded areas 
show th e uncertainty at T = 4 x lO^K of the esti- 
mate of iFisker et"al1 (|2006f ) and lTan et all (|2007[ ) 
respectively. 



3.1. High values of / 

For log(/) > —1.5 burning is unstable for ac- 
cretion rates up to log(/) =0.2-0.4 (see Fig. [T]) 
and a H/He triggered burst occurs after the 
i50(a,7) ^^Neif3+,v) ^^F{p,a) ^^0{p,j) i^F(p,7) 
'^^Ne(0+[v) ^^F(iP a)^'^0 cycle (for details see 
Cooper fc Naravanl[2006al : lFisker et al.ll2007t ) has 
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been active for several hundreds of seconds before 
the runaway. Once the 3a reaction triggers, the 
-'^^Ne(p, 7) ^°Na breakout reaction starts the re- 
process. This activates the ^^0{a,p) ^''F-reaction 
and creates a direct reaction flow from hehum into 
the rp-process viz. from 3a to ^^C(p, 7) ^'^N(p, 7) 
"0(a,p) ^'^F(p,7) iSNe(a,p) ^iNa. This reaction 
circumvents the T1/2 = 76.4 s /3+-decay of ^"'O 
and allows a flow directly into the rp-process. 

For log(/) > —1.0 burning becomes stable for 
log(/) =0.2-0.4 independent of the exact value of 
the ^^0(a,7) ^^Ne reaction. This means that pre- 
vious theoretical and com putational estimates of 
the critical accretion rate ("Fushiki fc Lambl llQSTt 



Bilds tcn IQQSt iFisker eTal. 2003; Wooslev et al 



iooi iHeeer et al.ll2005l) which were based on the 



rate of lCaughlan fc Fowleil (|l988l ). i.e. log(/) = 0, 
are independent of the exact value of the reac- 
tion rate. For log/ —1.5, small uncertainties 
the -'^^0(a,7) ^^Ne reaction rate would be suffi- 
cient to change the critical a ccretion rate by al- 



most an order of rnagnit ude (jCooper &: Naravan 



2006al : iFisker et al.ll200d ). 



3.2. Low values of / at low accretion rates 

For log(/) < -2 and log(/) < -0.5, the 
"'^^0(a,7) ^^Ne-reaction is too weak, so the above 
mentioned cycle never activates the ^*0(a,p) 
-'^'^F-reaction; the small amount of unstable he- 
lium fizzles, and the thermonuclear runaway never 
happen s. This res ults in an oscillating behavior 
fFisker et al.ll200d ). 



Cooper fc Naravan ( 2006a[ l speculated that a 



low value of / could stabilize bursts at lower ac- 
cretion rates. Using the time-dependent model at 
low values of /, wc find the same stable/unstable 
transition accretion rate for the helium trigger at 
lo g(/) ~ —0.5 thus comphmen ting the calculations 
of Cooper fc Naravai] ( 2006a ). However, It should 
be noted that for log(/) < —2 the instability of the 
helium trigger never results in a burst as explained 
in the preceding paragraph and the next subsec- 
tion. Therefore this behavior is marked with an s 
in Fig. m 

3.3. Low values of / at high accretion rates 

For log(/) < -2 and log(/) > -0.5, the 
higher accretion rate increases compressional heat- 
ing which heats the atmosphere to several hundred 



million degrees. Despite the low value of / the 
natural temperature dependence of the ^^0(a,7) 
^^Ne-reaction means that helium burning triggers 
the rp-process. However, since the value of / is 
small, the rp-process is fed very slowly resulting 
in increasingly longer luminosity rise times on the 
order tens of seconds as / decreases with burning 
becoming stable at very low values of / as shown in 
Fig. [21 For these values of /, the accretion rate at 
which unstable burning transitions to stable burn- 
ing is still uncertain. 

Fig. [2] shows luminosity peaks with an almost 
symmetrical decay suggesting that the atmosphere 
is driven by an oscillating burning front rather 
than a thermonuclear burst followed by a ther- 
mal decay. Such long rise times have never been 
observed at any accretion rate, so this result 
strengthens the argument for a lower limit on 
the ^^0(a,7)^®Ne rate but does not change the 
bound s on the astr qphysically based lower limit 
set by IFisker et al.l |2006l . This means that the 



oscilla tory burning behavior found in IFisker et al 
( 2006h is most likely not the generator of the 
^fC that fuels superbursts (|Cumming fc Bildsten 
200 1| ). although such a generator has not yet 
been found in other multi-zone simulations either 
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Fig. 2. — This figure shows the luminosity for 
three different values of / as a function of time 
for simulations with an accretion rate of log(J) = 
—0.4. Notice the very long rise times and the sym- 
metrical decays. 
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3.4. Comparison to observations 

Observations indicate that type I X-ray bursts 
occur for ]ne(T] < — D.fi and cease for ]ne(T) > 



-0.5 (Ivan Paradiis et al.lll979. 


1988; Cornelisse et al. 


2003;lRemillard et al.ll2006l; 


Gallowav et al.ll2006l). 



In contrast with the global line ar stability analysis 
of ICooper fc NaravanI (|2006al) , our full network 
time-dependent simulations show no single trial 
value of / that is consistent with observations over 
the entire range of accretion rates we consider. 
This suggests that more complex physics must be 
included in the simulation. The most important 
parameters of the bursting behavior as a function 
of accretion rate are the accreted composition - 
here assumed to be solar - and the mass and radius 
of the neutron star. Other possible factors include 
the sedimentation of CNO ions dPeng et al.li2r)07h 
as well as the geometry of the flow from the accre- 
tion disk thro ugh the boundary layer and o nto the 
neutron star ( Inogamov fc Sunveav 1999t ). Here 
the flow of the accreted matter may change as the 
thermonuclear burning transitions from unstable 
bursts to stable burning. 

4. The new ^^0(a,7)^^Ne reaction rate 

The ^^0(a,7)^^Ne reaction rate is dominated 
by the resona nce contributions at t emperatures 
above 0.1 GK ( Langanke et al. 19861 ). However, 
the partial widths Fq and of the relevant states 
with excitation energies of 4-5 MeV are not suf- 
ficiently well known for calculating the reaction 
rate reliably. Of particular interest is the res- 
onance at an excitation energy of 4.03 MeV in 
"'^^Ne, just above the ^^0-l-a threshold. This res- 
onance dominates the reaction rate at tem pera- 
tures below 0.6 GK ( Langanke et al. 19861 ) and 
therefore determines the temperature conditions 
for the breakout from the h ot CNO cycles. Its 7 
width was first mea sured by Tan et al. J 2005[ ) and 
later confirmed by iKanungo et al .1 (|200()h . How- 
ever, its s mall a-decay branching ratio of the or der 
lO"'' (;see lDavids et al.ll2003l:[Rehm et al.ll2003l for 
previous upper limits) prevented a reliable esti- 
mate of the reac tion rate until the recent work by 
Tan et al.l (|200<t . 

A new ^^0(a ,7)^^Ne rate was proposed by 
Tan et al. I (I2OO7I) based on the measured a-decay 
branching ratios and lifetimes of the relevan t 
states in ^^Ne (for details, see iTan et all l2007t ). 



This rate is shown in Fig. |3| with one sigma un- 
certainty indicated by the darker area, and it is 
significantly improved as compared to the pre- 
vious model-dependent estimate by iFisker et al. 
(I2OO6I) . In addition, dominant contributions at 



higher temperatures from the states at 4.14 and 
4.60 MeV were unexpected. This new rate not 
only allows a better identification of the ignition 
conditions of X-ray bursts but also permits the im- 
proved analysis of the dynamics and mechanism 
of X-ray bursts as demonstrated below. 
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Fig. 3. — Figure shows the newly measured reac- 
tion rate along with la upper and lower experi- 
mental u ncertainties. The u ncertainty range dis- 
cussed by iFisker et al. ( 20061 ) is also shown. 



4.1. Astrophysical consequences 

Rep eating the calculations of IFisker et al 



(|2006f ) for the lower limit of the ^^0{a,'^y^'Ne 
reaction rate. Fig. |4] shows a comparison between 
the luminosity as a function of time for the previ- 
ous lower limit and the present experimental lower 
limit for a constant accretion rate. The results 
show that the new rate within its experimental 
uncertainties is sufficient to trigger the observed 
sequences of bursts for a constant accretion rate. 

By doing a parameter study covering the pre- 
vious experimental uncertainty range for the 
^^Ofg, 7)^^Ne react i on it was shown in section |3| 
and in IFisker et all (|2006t ) that past X-ray burst 
models have been subject to a very large uncer- 
tainty due to the uncertainty of the ^^0(a, 7)^^Ne 
re action rate. Fortunately, t he widely used rate 
of lCaughlan fc Fowled (Il988h (viz. log(/) = 0) is 
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Fig. 4. — This figure shows the luminosity as a 
function of time for a simulation with an accretion 
rate of M = lO^^gs^^ comparing the lower limit of 
the "'^^0(a, 7)^^Ne for this work vs. t he lower limit 
of the ^^0(a,7)^^Ne rate used by iFisker et al 



( 20061 ). We note that the reduced uncertainty of 
the rate presented in this work corroborates obser- 
vations and thus constitutes a major improvement 
of the rate. 



very c lose to the newly measured rate o f iTan et al 
POOTI ). This means that the conclusions of ear- 



lier c a lculations (iFushiki fc LambI 1987 Bildsten 



19981: iFi sker et al.l I2OO3I: IWooslev et al.' |2004 



Heger et al. 2005: Cooper fc Naray an 2006a) hold. 



Since the ^^0(a, 7)^^Ne-reaction rate now has 
experimentally determined upper as well as lower 
limits, we can determine the accuracy of the theo- 
retical estimate of the critical transition accretion 
rate between the steady state burning phase and 
the burst phase in accreting neutron star binary 
systems. 

We already noted from Fig. [1] (see subsection 
13. ip that log(/)crit. = 0.2-0.4. This cor roborates 
the rnajority of previou s simulations (iRembged 



I999I : IFisker et al.l 120031: iHeger et al. 2005) and 
calculations (jFujimoto et al .1 119811 : lBildsten.l998h 



which have determined the transition point to be 
around M ^ 2.1 x lO^^gs"^ for a fiducial neutron 
star with i? = 10 km and M = 1.4M0 and adopt- 
ing a solar composition for the accreted matter. 



Using the rate of lTan et al.l (|2007l ). several sim- 
ulations were run for different accretion rates in 
the log(/)crit. =0.2-0.4 range while tracking the 
luminosity resulting from the nuclear burning. 



The results shown in Fig. [5] show that the burning 
becomes stable for M > 1.9 x lO^^gs"^. Identical 
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Fig. 5. — We used the newly measured rate to 
calculate the luminosity originating from the nu- 
clear burning as a function of time for different 
accretion rates. As is seen from the constant lu- 
minosity on the graph, the burning is stable for 
M > 1.9 X lO^^gs-^ 

simulations were performed using the one sigma 
upper and lower limits of the newly measured re- 
action rate as shown in Fig. [31 While the up- 
per limit yields the same transition accretion rate, 
the lower limit increases the transition point to 
M « 2.1 X lO^^gs"^. The astrophysical uncer- 
tainty in the determination of the accretion rate 
at the transition point has thus been reduced to 
less than 10% compared to previous rate-induced 
uncertainties of one order of magnitude in the ac- 
cretion rate c.f. Fig. [T] and Fisker et al. ( 2006f) . 



Conclusion 



Three important points: 



1. We corroborate the stab ility analysis of 



Cooper fc NaravanI (|2006al ) showing that the 



atmosphere is stable towards runaways if the 
"'^^0(a,7) ^^Ne reaction rate is low. How- 
ever, we also show that this instability does 
not lead to observable bursts if the ^^0(a, 7) 
^^Ne reaction rate is low. 

2. The new measurement is close to the previ- 
ous an d widely used rate of lCaughlan &: Fowler 
(Il988l) . so the conclusion s of previous X- 



ray burst simulations (e.g. Fushiki fc LambI 
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1987t iBildsteiil Il998l: IPisker et all l2003t 



Wooslev et alfeOOilHeger et al.l2005tlFisker et al 
do not change. 



3. The new measurement of the ^^0(a,7)^^Ne 
reaction rate significantly reduces the model 
uncertainty but does not result in a value 
that is in accord with astronomical obser- 
vations. Therefore further studies of the 
other determinants (mass, radius, accretion 
composition, neutron star core, sedimenta- 
tion, and possibly the accretion geometry) 
are needed. Such studies are currently un- 
derway. 

The simulations in this paper demonstrate why 
it is important to consider the uncertainty associ- 
ated with the input parameters of an X-ray burst 
simulation as it can significantly influence the pre- 
dicted observables. Furthermore, they show how 
experimental nuclear data can complement obser- 
vational results for a better understanding of the 
complex interplay between the fuel supply and 
burning processes on the surface of accreting neu- 
tron stars. 
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discussions and the referee for helpful suggestions 
that improved the clarity of the paper. This work 
is supported by the National Science Foundation 
under grant No. PHYOl-40324 and the Joint In- 
stitute for Nuclear Astrophysiclll, NSF-PFC under 
grant No. PHY02-16783. R. L. C. is supported 
in part by the National Science Foundation under 
Grant No. PHY99-07949. 
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